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We investigate the interplay between the formation of lattice and magnetic polaron in the case 
of a single hole in the antiferromagnetic background. We present an exact analytical solution of 
C ' the Holstein-t-J model in infinite dimensions. Ground state energy, electron-lattice correlation 

function, spin bag dimension as well as spectral properties are calculated. The magnetic and hole- 
lattice correlations sustain each other, i.e. the presence of antiferromagnetic correlations favors the 
formation of the lattice polaron at lower value of the electron-phonon coupling while the polaronic 
effect contributes to reduce the number of spin defects in the antiferromagnetic background. The 
crossover towards a spin-lattice small polaron region of the phase diagram becomes a discontinuous 
transition in the adiabatic limit. 
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I. INTRODUCTION 



. . . . 

^ ' A single hole in an antiferromagnetic background represej^s a widely studied problem in solid state physics due to 
'"^ , its relevance to the problem of high- Tc superconductivity.Em High- Tc compounds indeed show an antiferromagnetic 
' undoped phase which is gradually degraded and finally destroyeidJpy doping. The electronic properties at low doping 
O ' are therefore often described in terms of Hubbard or t-J modelJj'B 

The problem of the motion of a single hole, far from being a pure academic issue, can be representative of the 
extreme low doping case of these materials. It is worth to note that even in the case of a single hole the problem of 
I intermediate/strong magnetic interaction with an antiferromagnetic background is a non trivial many-body problem. 
K*" ' The difficulty consists in describing the dressing of the hole by a. doud of spin background excitations (the magnetic 
or spin polaron) which can coherently moves as a quasiparticleM~tl The situation is similar to that of a small lattice 
■ polaron i.e. the case of an electron moving together with a phonon cloud which represents the lattice deformation 
induced by the presence of the charge. 3 

The connection between lattice and magnetic polaron goes however beyond a methodological interest. Thercj-are 
fSl indeed several observalipns of a sizable interplay between electron-phonon and magnetic interaction in cupratesO as 
. well as in manganites.Ed Purpose of this paper is to explore in detail the physical consequences of this interplay, in 
' particular with regards to the lattice and magnetic polaron properties. To this aim a non perturbative way is clearly 
^ , needed. 

^ ' We present an exact analytical solution of the Holstein-t- J model for a single hole in infinite dimension. Ground 
state energy, electron-lattice correlation function, spin bag dimension as well as spectral properties are calculated. 
^ ■ We find that the lattice (spin) polaron formation depends strongly on the magnetic (hole-phonon) interaction. We 
Q identify thus regions of phonon assisted magnetic polaron as well as magnetic induced lattice polaron. The extension 
^ of these regions are strongly dependent on the adiabatic ratio and they vanish in the adiabatic limit. In that regime 
lattice and magnetic polaron formations are strictly tied each other. These general results could help to explain the 
strong interplay between lattice and spin degrees of freedom in cuprates and in manganites. Finally we discuss the 
major drawback of our approximation and we give some ideas to overcome it. 



II. THE HOLSTEIN-t-J MODEL AND ITS DMFT SOLUTION 

Let us consider the Holstein-i- J model defined by the HamiltonianiElIll 
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where c]^ are the electron operators in the presence of infinite on-site repulsion that prevents double occupancy 
[cj^ = c|^(l — n^a)], the phonon operators and S^''^'~ respectively the z component and the raising and lowering 
spin operators. The first term in Eq. (|l]) describes the hopping of the electrons on nearest neighbors of a square 
lattice, the second one the direct and the exchange interaction, the third one the local electron-phonon interaction 
coupled to charge density, and the fourth the Einstein phonon frequency. The choice of a hopping matrix element 
equal to i/2 gives rise to a band with bare bandwidth t. The model can be straightforward generalized in infinite 
dimensions by using the usual rescaling: t — > t/^/z, J J/z, where where z is the coordination number. For a 
hypercubic lattice the coordination number is 2: = 2d, while for a Bethe lattice z = d. 

All through out this paper we shall consider one hole created on the anti-ferromagnetic half filled state. Antifer- 
romagnetic state is described in terms of a classical Neel ground state. A convenient approach to this aim is the 
spin wave theory applied to lattice model as can be found e.g. in Ref. A useful effective Hamiltonian can be 



thus derives by mainly following the discussion in Ref. 11, generalized now in the presence of a Holstein electron- 



phonon interaction. The Hamiltonian is first transformed by a canonical transformation into a ferromagnetic one. 
Then "hole" and "spin" defect operators are introduced, respectively as fermionic h and bosonic a operators on the 
antiferromagnetic ground state. The resulting Hamiltonian reads thus: 
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where we have neglected the hole-hole terms since we are interested in a single hole in an antiferromagnetic background. 

Eq. (|^) can be significantly simplified in infinite dimension. In that limit indeed the two terms of the fourth line can 
be shown to be negligible since they contribute only at 0{l/d). In addition, in the absence of any boson condensate 
(a), (a^), which should destroy the antiferromagnetic background and which is forbidden in our context, also the last 
two terms of the third line can be dropped. We end up thus with the effective Hamiltonian of the Holstein-f- J model 
valid in infinite dimension: 

-gJ2hlh,ib. + bl)+LOoY.blb, 

i i 

The first term of Eq. (^) describe the kinetic hopping of one hole on the antiferromagnetic background, which is 
accompanied by the creation (destruction) of a spin defect which breaks (restores) 2z magnetic bonds with individual 
energy J/4z. In addition we have the usual local electron-phonon interaction which couples the hole density to the 
local phonon. The last term in Eq. (||) can be absorbed in the definition of the hole chemical potential which, for the 
single hole case here considered, has to be set at the bottom of the hole band. It is important to note that, although 
written in terms of hole operators, the phonon part (free phonon part -I- hole-phonon interaction) is formally identical 
to the Holstein model. This is not a trivial result since in a half-filling case all the electrons are coupled with the 
phonons and one should in principle deal with a many-body problem. As a consequence we are thus able to reduce 
the many-body problem to a single-particle (one hole) system interacting with phonons and with spin defects through 

Eq. ®- ^ ^ 

Although the Hamiltonian (H) looks much more affordable than (|^), the analytic study of its properties is still a 
quite hard task at finite dimension. The problem can be much simplificd_h,owever in infinite dimensions where, as 
we are going to see, the limit of infinite coordination nujaber z — oo,tj1l3 all together with the retraceable path 
constraint enforced by the antiferromagnetic background ,113 provides an exact solution. A explicit derivation can be 
found in Appendix In this section we only summarize the final equations which determine in a self-consistent way 
the hole Green's function. 
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A crucial point is the possibility of writing the self-energy of the local propagator as the sum of two contributions, 
labelled as Shop(w) and Soi_ph(w), which closely resemble the functional expressions of the hopping and phonon 
self-energy respectively in the pure t-J and Holstein models, but which are now evaluated in the presence of both 
exchange and phonon interactions. We can write thus: 

^(^) = — ^ I I — r-v (4) 



where the "hopping" contribution is given by0 



Shop(c^) = jG{uj - J/2), (5) 



and the "phonon" self-energy can be expressed by means of a continued fraction 

G^^{uj - uja) - 



Eel-ph(w) = , (6) 
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Gf. ^{uj — 3uJo) — . . . 



where 



Gr'H-^- J/2). (7) 

It should be stressed again that both I]hop(^) and Soi--ph(w) arc functions of the total Green's function G [Eqs. 
(^"(0)] which contains the full dynamics (hopping, exchange, electron-phonon) of the system. This not trivial self- 
consistency accounts thus for the complex interplay between the phonon and spin degrees of freedom. 

Eqs. (U)-(|7|) represent a closed self-consistent system which we can be numerically solved by iterations to obtain 
the explicit exact expression of the local Green's function G{uj), and hence any local one-particle relevant property of 
the system. 

The formal scheme looks quite similar to the dynamical iaean field theory in infinite dimension for a Bethe lattice, 
applied for instance at the purely electron-phonon system.tS However, due to the antiferromagnetic background, the 
physical interpretation is quite different. 

Due to the orthogonality of the initial and final antiferromagnetic background, the noUpLpcal component of the 
Green's function in the Holstein-t-J model for J = is strictly zero Gij{uj) = G{uj)dij jtli whereas for the pure 
Holstein model Gi^j{uj) is finite and provides information about the non-local dynamics: G{\i,uj) ~ l/[ct; — ek — S(w)]. 

In addition, the magnetic ordering has important consequences also on the local Green's function Gii{Lu). In the pure 
Holstein model for instance Gu^lo) takes into account any generic dynamics which occurs back and forth a given site 
whereas in the Holstein-t- J model the electron must follow a retraceable path in order to restore the antiferromagnetic 
backgroundO A Bethe- like dynamics is thus enforced by the magnetic ordering regardless the actual real space lattice. 
The object made up by the hole plus the local modification of the spin configuration due to the presence of the hole 
is the "spin polaron" . 

The local constraint Gij{uj) = G{Lu)Sij induced at d = oo in the Holstein-t-J model by the antiferromagnetic 
background can appear a quite strong simplification. However, it should be reminded that it holds true as long as the 
antiferromagnetic spin configuration can be assumed to be frozen, in particular, as long as the Hamiltonian does not 
induce spin dynamics. This-|is the case of the effective Holstein-t-J Hamiltonian in infinite dimension where spin 
fiuctuations are neglectedO The existence of the spin polaron itself could be questioned at finite dimension where 
spin fluctuations are operative. However, several numerical and analytic studies have shown that the restocing of spin 
fiuctuations does not destroy the spin polaron object, but opens coherent channels of hole propagation. B"I3 In this 
situation the spin polaron can thus propagate as a whole through the crystal. 

Many studies have investigarUed, the motion of the spin polaron and determined its k-|disp|epion and optical con- 



ductivity, both . in the absenceOO EI and in the presence of electron-phonon interaction. t3'Eirc3 However, apart few 
exceptions ,tjO the internal degrees of freedom of the spin (lattice) polaron object have not so far been much in- 
vestigated. In the present work we mainly focus on the formation of the spin polaron and on its internal structure 
properties (size, binding energy etc. . . ). These quantities are local features which, we believe, are only weakly affected 
by the itinerant nature of the spin polaron. In this perspective we think that the infinite dimension approach here 
considered provides valuable information on the spin polaron formation and on its interplay with the local (Holstein) 
electron-phonon interaction. 
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III. SPECTRAL PROPERTIES 



In order to have a complete description ofJie physical properties of the Holstein-t- J model, it is useful to identify 
three independent dimensionless parameters :tZI the adiabatic ratio wo/i, the electron-phonon coupling A — jujot and 
the exchange interaction J jt. Another important parameter to be defined is the multiphonon constant a = g/ojQ that 
is a local quantity which does not involves electron hopping t. Limiting case are A = where the system reduces to 
the t-J model, and J /t = where the Holstein model on an antiferromagnetic background is recovered. 

A naive look could regard this problem as a simple interplay between two energy scales: the exchange J/t which rules 
the magnetic properties, and A related to the electron-phonon coupling. The dominance of one of them would therefore 
determines the overall properties of the system, while the weaker one could be considered as a perturbation. However, 
as we are going to see, this picture is too simplistic. A more accurate description of the physics must take into accpunt 
first of all the role of the adiabatic ratio (ujo/t) which in the pure Holstein model rules also the polaron crossover .c3 In 
particular, we can expect that the magnetic (lattice) polaron formation induces a drastic renormalization of the kinetic 
energy scale t. The "effective" adiabatic ratio will depend in an implicit way on the electron-phonon and magnetic 
interactions. Then the magnetic (lattice) polaron formation induces a drastic renormalization of the kinetic energy 
scale which in its turn affects the phonon (magnetic) properties and eventually leads toward intermediate/strong 
hole-phonon couplings. It is interesting to note that the multiphonon parameter a = g/ujQ does not depend on t, so 
that it can be considered unaffected by local hopping renormalization (we remind that phonon frequency coq is not 
screened by electron-phonon interaction for the single hole case). 

Let us now first discuss spectral properties by studying the spectral density defined as A{u!) = — (l/7r)Im[G(ijj)] 
directly accessible by the knowledge of Green function. 

All through this paper we consider a Bethe lattice a the standard semicircular density of states with bandwidth 2t. 
It should be noted however that, since the antiferromagnetic background enforces a retraceable path approximation, 
a semicircular density of states is recovered independently of the chosen crystal lattice, i.e. also for a hypercubic 
lattice. In other words Eqs. are valid for any lattice structure provided the coordination number is infinite. The 
assumption of a non-retraceable Bethe lattice however allows to classify explicitly all the hopping processes which 
leads to the magnetic polaron formations Jill Moreovepin the Bethe lattice, contrary to hypercubic infinite bandwidth 
case, the single polaron problem can be well defined.ll3 To obtain the spectral function we have iterated Eqs. (0)-(0) 
up to numerical convergence by using a truncation in the continued fraction Eq. (|^) according to the procedure 
outlined in Ref. |l^ i.e. we truncate the continued fraction Eq. (^ at a-stage iVph ^ a^. Moreover we notice that a 
continued fraction arises also in the absence of hole-phonon interactionO once Gt (Eq. (§)) is substituted in Eq. (Q). 
We have found that in order to properly generate a certain number M of the magnetic poles of the spectrum we have 
to choose a truncation Nj ^ M. 

In Fig. |l| we show the evolution of the spectral density for moderate values of J/t and A — 0.5, LUo/t = 0.5. For 
J/t = the spectral density is made by a continuum with incipient structures due to the electron-phonon coupling. 
This shape is characteristic of-weak electron-phonon coupling in a small-intermediate adiabatic regime oj^/t 1, with 
no well defined polaron peakJl3 By switching on the exchange interaction J/t = 0.2 the continuum spectrumjs split 
in a set of magnetic peaks. This trend is quite similar to what happens in the infinite dimensional t-J modelcZI with 
an additional modulation due to the underlying electron-phonon features. We can thus think of the resulting spectral 
function as ruled by different couplings on different scale, where the gross structure on scale t is determined by the 
pure electron-phonon interaction superimposed by the fine structure on scale J given by the magnetic peaks spaced as 
(J/t)2/3 (for small J/t). By further increasing J/t we pass an intermediate regime where magnetic and phonon peaks 
are mixed together and eventually we have a purely phononic spectrum with a set of peaks starting from — g^/cja-and 
equally spaced by ujq- This is indeed just the characteristic spectrum of the atomic Holstein model with — 0.5.E3 As 
discussed above, the origin of such a behavior is the strong hopping amplitude renormalization due to the magnetic 
polaron formation for J > 1, which drives the system from an almost adiabatic case to an effective antiadiabatic 
regime LOo/t* 3> 1. It is however surprising that a strong magnetic interaction J/t ^ A yielded a purely phononic 
spectrum!. 

In order to better understand the global evolution of the spectral function it is interesting to look at what happens on 
a larger energy scale. In Fig. ||we have thus plotted the the behavior of A{u!) on a linear-log scale for J/t = 0, 5, 10, 20. 
We see that the spectrum for finite J/ 1 is roughly determined by replicas of the J/t = spectrum equally spaced by J 
and with vanishing spectral weight. On a closer look we find that the structure can be derived from the spectrum of 
the pure t-J model, made for J/t ^ 1 by equally spaced magnetic peaks ^„ a„J((jj — nJ/2), by broadening each peak 
according the electron-phonon interaction as in J/t = case. The inverse occurs in for J/t < A. In this case we found 
magnetic structures on a small energy scale J coexisting with a phonon structure spread on energy t. On the contrary 
the spectral function for large J/t can be thought as built by pure magnetic structure for large energy superimposed 
to a finer phononic structure made by peaks separated by the bare phonon frequency ojq. This feature resembled the 
"interband" transitions found in Ref. ^ in the anti-adiabatic regime. Of course when the magnetic coupling strictly 
goes to infinite J/t oo the high energy magnetic peaks are shifted to infinite energy and completely loose their 
spectral weight. The spectral function reduces in this case to what shown in the lower panel of Fig. |l|. 
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IV. GROUND STATE PROPERTIES 



Ground state properties can be derived in a direct way by the knowledge of the Green's function which aUow the 
evaluation of the ground state energy Eq through the determination of the lowest band edge or the lowest pole. 

The character of the ground state can be deteriBiped by the knowledge of several relevant quantities that can be 
evaluated using the Hellman-Feynman theorem as:E2l 

i) The mean number of phonons A^ph = A^ph = dEo/dajQ; 

ii) The mean hole-phonon correlation function Cq — {h^h{h + 6^)): Co = —dEo/dg; 

iii) The mean number of spin defects Ng.d. — (a^a): N^.d. — 2dEo/dJ; 

iv) The effective hopping amplitude t*: t* /t = 2d\Eo\/dtM 

Quantities i) and ii) shed light on the lattice polaron formation process. A sharp increase of Cq {Np\^) is expected 
around some intermediate value of the hole-phonon coupling A in the adiabatic regime. Cq increase from zero to its 
strong coupling limit 2a-La?). The transition becomes a crossover which becomes smoother and smoother approaching 
the antiadiabatic case.liSEa ■— ■ 

Quantity iii) provides information on the size of the magnetic polaron.Ea In fact, since the retraceable path ap- 
proximation is enforced by the antifeptHiaagnetic background in infinite dimension, it is clear that the A's.d. gives 
the length of the string of spin defectsE30tZl which, in the Bethe lattice, is also the size of the magnetic polaron. In 
the zero exchange limit J — )■ no energy cost is associated with a spin defect and the size of the magnetic polaron 
diverges (large magnetic polaron limit). In the J — > oo case instead spin defects are unfavored and the magnetic 
polaron becomes pretty local (small magnetic polaron limit). It could be appear surprising that information of non 
local quantities, as the magnetic polaron size, could be available in the local approach we are using. However it should 
be reminded that this is only a mean quantity. The price we are paying by using the local approximation exact in 
infinite dimensions is the impossibility to have information on the probability distribution to find a path with a given 
length n. We shall see later that a careful inspection of the dependence of N^a. on the magnetic coupling J can 
nevertheless provide valuable information about the rough shape of the magnetic polaron size distribution. 

In order to establish however a criterion for the large-small magnetic polaron formation, we define the value 
-^s.d. = 0.5 as the conventional transition between large and small magnetic polaron. We can indeed think that 
for Ns.d. > 0.5 the probability to find paths with length n > 1 is larger that the probability to have a completely 
magnetically trapped hole (n = 0), and we are therefore dealing with "large" magnetic polarons. 

Finally the tendency to localization that can be due either to magnetic or hole-phonon interaction can be deduced 
from the behavior of the effective hopping iv) . 

Returning to the analysis of the lattice nolfljpjxxaae we could certainly use as well the criterion A^ph > 1 or Cq > 1 
to identify the lattice polaron formation.li3£J't§EZlc3 However, because of the local nature of the electro n- rah onon 
interaction, we can gain a deeper insight by looking at the probability distribution of the phonon numbers .Eilc3 The 
probability distribution of phonon numbers is defined as 

P{n) = I {n\h\0) p (8) 

where |n) is a state with zero holes and n phonons and |0) is th e gr ound state of the single hole. It can be obtained as 
a residue at the ground state enecgy uj ^ Eq of G""(cij) [Eq. ( |A3| )]. From the same procedure outlined in Appendix 
^ (see also Ref. ^6|) one obtainsH 

G™ (c.) = \ ^ , , ^ , , , (9) 



where Y^^.^(lo) represents the processes related to the emission of a phonon from the state |n) 

+ 2)g 



Scm(w) = , „^ 2 ' (10) 



^{uj — ncjQ ~ ujq) — 



Gf. "'"(w — nojQ — 2ljq) 



{n + 3)g2 



G( ^{uj — nujQ — Swo) 



and Sabs('^) takes into account the absorption processes which are allowed also at zero temperature by the initial 
71-phonon state 
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The propagator G^ ^{uj) is defined in Eq. (Q), so that G""(cj) is a direct by-product of the self-consistent solution of 
Eqs. From that we immediately obtain the phonon number distribution as 

P,„,.f&riy , ,,2) 



For a noninteracting system (17 — 0) the phonon distribution contains only a (5-peak at n = 0: P{n) — <5n.o. 
Switching on the electron-phonon interaction, the onset of local lattice distortions are reflected in a shift of total 
weight towards higher multiphonon peaks. We can now unambiguously identify the lattice polaron formation with 
the condition 

P{n = 0)<P{n=l). (13) 

It is easy to check that this definition reproduces the well known results for the Hjalsicin model, namely the criterion 
A > Ac in the adiabatic limit (Ac — 0.844), and a > 1 in the antiadiabatic one.EZrCjll3 It could be worth to remark 
that the polaron transition occurs however in a different way in the two limits. In the nonadiabatic regime the most 
probable phonon number n evolves in a smooth way by increasing A from n = to higher numbers n ~ a^. The 
dependence on A becomes sharper and sharper by decreasing the adiabatic ratio ujo/t and in the adiabatic limit 
uo/t = n jumps in a discontinuous way from h = for X < Xc to n = 00 for for A > Ac. 

In Fig. ^ we plot the transition curves corresponding respectively to P{n = 1) > P{n — 0), P{n = 2) > P{n — 1), 
P{n = 3) > P{n — 2), etc. in the X-J/t phase diagram for ojQ/t = 0.5. Lattice polaron formation occurs on the left 
line corresponding to P{n = 1) = P{n = 0). We notice a strong dependence of the lattice polaron formation on the 
magnetic energy. In particular increasing the exchange coupling J/t the lattice polaron formation is shifted to smaller 
values of A until at J/t — > 00 lattice polaron formation is ruled by the antiadiabatic criterion = 1. 

The role of the magnetic interaction in driving the system towards an effective antiadiabatic limit is even more 
evident when we draw the lattice polaron phase diagram in the X-uJo/t space (Fig. ^). For zero exchange coupling 
J/t = 0.0 (solid line) it is possible to distinguish an adiabatic regime, where the lattice polaron formation is ruled by 
the condition A ^ 1, and an antiadiabatic regime where the polaron occurs for > 1 By switching on the magnetic 
interaction the polaron crossover approaches the line = \ and the validity of the antiadiabatic criterion is extended 
for smaller values of Lo^/t. 

We can now address the open issue concerning the modification of the lattice polaron criterion in the presence of 
electronic correlation. The point is to determine whether the simple relation A < Ac in the adiabatic regime could 
be generalized by introducing properly scaled parameters. Two alternative pictures have been debated in literature. 
According the first onctj the relevant parameter in the presence of electron-electron and magnetic interaction is 
the ratio between the lattice polaron energy g^/wo and the purely electronic ground state energy in the absence of 
hole-phonon interaction -Eo.mg = -£'o(A = 0) 

K = IP I - (14) 

An alternative point of vie w0 regards the effective hopping amplitude t*^^ — t*{X = 0) as the main renormalization 
effect of the exchange coupling 

According these two ideas the relation A = Ac should be replaced in the presence of magnetic interaction by A* = Ac. 

We have carefully checked the validity of these two criteria within our exact solution in infinite dimension. We 
found that both of them fail to locate correctly the polaron crossover in the presence of magnetic coupling because 
the effective adiabatic ratio is increased by the decrease of kinetic energy (tjjjg) due to magnetic localization. This 
drives the system toward an anti-adiabatic regime in which the polaron crossover is ruled by the multiphonon constant 
a which is not renormalized by magnetic coupling. A naive way to take into account the reduction of the kinetic 
energy is to renormalize also the adiabatic ratio uja/t in similar way with Eqs. ([l^)-(|l5|), respectively a;o/|i?o,mg| 
and wo/l^mgl- We have checked also this criteria and we found that the lattice polaron crossover in the presence of 
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magnetic interaction can not be described in a satisfactory way even within this scheme, although it provides a better 
agreement than the simple renormalization of A. This just means that the magnetic and lattice degrees of freedom can 
not be separated, namely that the single hole Holstein-t- J here considered can not be mapped in a simple Holstein 
model with renormalized parameters. 

After having analyzed the process of lattice polaron formation, we can now investigate the properties of the spin 
polaron in the presence of hole-phonon interaction. We have just seen that the hole-phonon and magnetic interactions 
are not in competition. On the contrary the exchange coupling results to favor the lattice polaron formation. On the 
same foot we can expect that similar arguments hold for magnetic polaron formation, namely that the electron-phonon 
trapping favors small spin polarons. 

In Fig. ^ we show the mean number of spin defects A's.d. as function of the exchange coupling J/t for various 
electron-phonon couplings A and different adiabatic ratios uj^/t. In the pure t-J model (solid line) the small-large 
magnetic polaron crossover is denoted by a change of slope at about J/t = 1 which separates a iVs.d. oc J^^ from a 
A's.d. C!C J~^/^ regime.c3 

This trend is qualitatively unaffected in the anti-adiabatic regime uo/t = 2.0 for the electron-phonon coupling 
here considered A < 2. In this situation the lattice polaron formation occurs as a smooth crossover with negligible 
localization. Large-small spin polaron formation is driven therefore only by the magnetic interaction without any 
significant interplay between lattice and spin degrees of freedom. We can schematize this scenario as a two phases 
transition, large magnetic polaron/small magnetic polaron. 

The scenario changes approaching the adiabatic limit. In the intermediate regime l/Jq/I = 0.5 we can distinguish a 
weak electron-phonon coupling regime (A < 1) where the magnetic interaction is still the only relevant energy scale 
for the large-small spin polaron transition, and the strong coupling regime (A > 1) where lattice polaron formation 
interferes with the magnetic one. The case A = 2 (dot-dashed line) is representative of this regime. For very small 
J/t ^ we recover the usual N^.d. — cJ"^!'^ behavior, characteristic of the large magnetic polaron. Effective 
electron-phonon coupling is not sufficient to give the lattice polaron localization and it solely gives a reduction of 
the prefactor c due to the renormalization of the local hopping amplitude. For larger J, 0.05 J/t J; 5, we find a 
remarkable decrease of the mean number of spin defects Ns.d.- The origin of such a decrease is not magnetic since in 
this region A's.d. depends only weakly on the exchange coupling J. We can identify this regime as a small magnetic 
polaron induced by lattice polaron trapping. Finally, for larger interaction J/t > 5 the magnetic energy becomes 
strong enough to overcome the lattice polaron localization and we recover a pure magnetic trapping. The overall 
evolution from large magnetic polaron to small magnetic polaron can be described as two crossovers: large magnetic 
polaron/small (magnetic) polaron induced by lattice polaron localization/small magnetic polaron. 

Crossovers become even more marked as approaching the adiabatic regime and become a discontinuous transition 
in the adiabatic limit {wo/t = 0.0 in Fig. |^). Detailed calculations for this particular limit (static lattice distortions) 
have been explicitly carried out in Appendix The dependence of A's.d. on J/t is drastically different in the weak 
A < Ac and in the strong A > Ac coupling cases. In this latter case in particular the particles are almost perfectly 
trapped the intermediate region of small magnetic polaron induced by lattice trapping extends towards J/t = 0. In 
the adiabatic limit, for A > 0.844 we have always small lattice/magnetic polaron (see Appendix 

In the above discussion particular care needs to be paid in distinguishing between the character (small/large) and 
the nature (magnetic of phononic) of polaron transition. We can find a similarity between the small/large magnetic 
polaron and the small/large lattice polaron transition: both of them depend on the local probability of the hole.ta Hon 
from a site i to an other site j. A suitable quantity to express this concept is the effective hopping amplitude i*.E30'LZl 
Small values of t* /t denote strong localization of the hole, regardless its specific origin (lattice polaron trapping or 
small magnetic polaron). 

In Fig. ^ we plot t* /t as function of the exchange coupling J/t (left panel) and of the electron-phonon coupling A 
(right panel) for the cases uJo/t = 0.0 and ujo/t = 0.5. Comparing the left panel of Fig. ^ with the corresponding cases 
in Fig. H we can clearly identify the trends discussed above. For A ^ 1 the electron-phonon interaction induces only 
a weak reduction of t*/t while the localization transition is essentially driven by the magnetic coupling J/t. When 
the electron-phonon coupling is strong enough (A 1) however the hole dynamics is strongly suppressed already at 
J/t = by lattice polaron trapping, and higher value of J/t are needed to further decrease the effective hopping 
amplitude by magnetic effects. 

This behavior is quite similar when we plot t* /t as function of A for different exchange coupling J/t (left panel). 
We note however that the decrease of the kinetic energy at uJo/t = 0.5 is steeper when induced by lattice polaron 
formation than by the magnetic one. This difference is amplified by approaching the adiabatic regime too/t <ti I. In 
the adiabatic limit wo/t = a discontinuous transition occurs at small value of J/t_around a critical value of the 
coupling which in the limit J/t = approaches the value found in the Holstein modelij (see Appendix 

We can now summarize the above study in a global polaronic phase diagram for the Holstein-t-J model, shown 
in Fig. ^ The solid line marks the lattice polaron formation according the criterion (13) and the dashed line the 
small/large spin polaron transition by A's.d. — 0.5. The dependence of the lattice polaron formation on the magnetic 
exchange J/t (solid line) and of the spin polaron transition on electron-phonon coupling A (dashed line) points out 
the strong interplay between the two kind of processes. In particular they are not competing but sustaining each 
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other. We can distinguish four regions characterized as foUows (see inset in Fig. |^): (A) no lattice polaron, large spin 
polaron; (B) lattice polaron, large spin polaron; (C) no lattice polaron, small spin polaron; (D) lattice polaron, small 
spin polaron. 

It is interesting to compare the evolution of the phase diagram with respect to the adiabatic ratio LOa/t. For large 
ujo/t the lattice polaron formation is not accompanied by a strong hole trapping but appear as a smooth crossover 
just like the magnetic transition. We can thus identify finite regions (B) and (C) where lattice and spin polaron can 
be established independently of each other. Approaching the adiabatic regime [uiq/I = 0.1) the phases (B) and (C) 
gradually shrink. 

Particular care is needed in the strict adiabatic limit ujo/t — 0. The criterion described Eq. ( p^ ) states the existence 
of a multiphonon state as a small polaron key feature. Due to the localized nature of the system, lattice distortions 
with vanishing quantum fluctuation are always present for any finite A (see Appendix ^) . This classical lattice state 
is indeed constituted by an infinite number of phonon giving Ac = by using the criterion of Eq. (p^. Nevertheless 
we can always identify a discontinuous transition from very small to large lattice distortions which survives up to 
J/t ~ 0.132. Explicit results are shown in Fig. ||. It is thus this transition which strongly reduces the electron 
local hopping t* and enforces the spin polaron. For larger J/t such a sharp transition disappears and also the small 
magnetic polaron formation becomes a smooth crossover (dashed line in Fig. ^). 

It should be stressed however that the extension of the strict adiabatic to realistic finite adiabatic ratio is relevant 
as far as quantum fluctuations are small with respect to the average lattice distortion. In the case of luq /t ~ 0.1 for 
instance quantum fluctuations are larger than the small average lattice distortion found in the adiabatic limit. In this 
case the multiphonon criterion as described by Eq. (|l3|) marks as well the small polaron crossover. 

V. CONCLUSIONS 

We have mapped the half-filled Holstein-t- J model on an antiferromagnetic background into a one-particle Hamil- 
tonian (hole interacting with phonons and spin defects) which can be exactly solved in infinite dimension in terms of 
a continued fraction. The method immediately gives access to the hole spectral density and ground state properties. 
The main results of our work can be summarized as: 

i) Magnetic and phononic excitations are well separated in the anti- adiabatic regime and/or in the strong magnetic 
regime. In this case phononic peaks are separated by the bare phonon frequency loq. 

ii) Magnetic and lattice correlations sustain each other. Not only polaron crossover is shifted to lower coupling by 
magnetic correlations as noticed in Refs. but also spin defects are reduced by polaronic effects. 

iii) Phonon retardation is affected for large J/t by magnetic correlations leading the system towards anti-adiabatic 
conditions in which the relevant electron-phonon coupling changes from A to a. Therefore it is not sufficient to scale 
A with the renormalized electron kinetic energy to locate the polaron crossover. 

iv) We identify a crossover region between regions of parameters space in which magnetic and lattice polaron 
occurs independently (B-C) and regions (A-D) in which they are mutually dependent. We find the rigorous adiabatic 
regime to be quite peculiar wherein the small magnetic polaron formation is always accompanied by a small polaron 
formation. 

The main drawback of our DMFT method lies in neglecting dispersion of spin waves which leads to hole coherent 
motion. This is evident in our spectra which are constituted by k-independent peaks. As a starting point to overcome 
this difficulty we explicitly included in Eq. (^ terms which would lead to hole delocalization as well as to spin wave 
dispersion in the next order in 1/z. A controlled way to include this processes is currently under investigation. We 
expect that taking into account the coherent quasiparticle motion of the hole due to the quantum spin fluctuation 
would modify the low energy features of the spectral function by giving rise to finite bands. However we think that 
the gross features of the spectral weight would not be strongly affected. In addition the existence of spin fiuctuations 
could lead to delocalization of the large polaron found in the adiabatic limit at weak coupling. 

We acknowledge useful discussion with M. Capone. 

The authors acknowledge support of the Italian Ministry of University and Scientific Research fund cofin-99. 

APPENDIX A: EXACT SOLUTION OF THE ONE HOLE HOLSTEIN-£-J MODEL IN INFINITE 

DIMENSION 

Let us consider the Hamiltonian in Eq. (^) which we write as H ~ Ht + Hl ■ Here Ht contains the hopping terms 
and all the other local contribution. In the absence of any hole, the ground state is just the antiferromagnetic one 
which can be written as \ AF) — |0; 0), where |0; 0) represents the antiferromagnetic background with no phonon and 
no spin defect on any site. In similar way we can introduce the notation 
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i(i),m(j),...;/,fc,...) 



In \/m 



a\al . . 



|0;0) 



(Al) 



to express the state with n phonon on the site i, m phonons on the site j , and defects of spin on the site Z, fc, 
Aim of investigation of this appendix will be the Green's functiontJ 



G,,M = (0;0 



1 



0;0 , 



which can be considered as the (0,0) element [G{uo) — G°"((jj)] of the generalized Green's function 

mi, 



In addition we introduce the atomic propagator 



which satisfies the following properties (i ^ j): 



0; rii, rrij 



0:ni 



hi 



uj-Hl 



TOi;0 ) , 



(A2) 



(A3) 



(A4) 



0; rrij 



1 



TT 

uj — Hl — najQ ■' 



rap I \ 

9i-i ~ n-uJo), 



(A5) 



and 



Si] Pi 



= 0:0 



1 



h, f h] 

h, ■ h\ 

^ UJ- J/2- Hl- poJo ^ 

g'^^{LO-J/2-pujo)dp^,. 



(A6) 

Eqs. (A^,A6) stem from the fact that the electron- phonon coupling in Hl is operative only on the site on which the 
hole stays. 

Let us now expand in Eq. (A2) the resolvent l/{uj — H) in powers of Ht: 

1 1 1 „ 1 



-Ht 



Lu — H Lu — Hl uj — Hl uj — Hl 
1 „ I „ I 



-Ht 



Ht- 



UJ — Hl uj — Hl uj — Hl 



(A7) 



It is easy to see that all the odd powers do not contribute in Eq. ( [A^ ) since they create or destroy a odd number of 
spin defects. More generally, since the initial and final states do not contain spin defects, it is clear that all the spin 
defects created in the dynamics must be destroyed on each site before to reach the final state. In infinite dimension this 
constraint selects only the retraceable paths. Each "step forward" is thus ruled by the term t / {2y^ J2 {a[3) f^a^ahp 
and each "step backward" by its complex conjugate. 

At the second order in t we have for instance: 



0; n. 



1 



-Ht 



1 



-Ht 



1 



{",/3> (7,<5) 



UJ — Hl uj — Hl uj — Hl 
1 t 



u-Hl2^ 



■mi] 
1 



co-Hl 



rrii', 



(A8) 
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Eq. ( |A^ ) can be expressed in terms of the atomic propagator g by introducing the identity operator in each hopping 
term: h^^aahp = ha'^a\M){M\hi3, hl^a\hs — '}2n h\\N){N\a\hs where |M), \N) are complete sets of states. It is 
easy to check that only the states = Hp \Pi\ Si){s^,Pi\ and X^at = Y.q Ift; Si)(si, %| give a non-zero 

contribution. Using the properties Eqs. (|A5|)-( A6) we have thus: 



(A9) 



Applying a similar procedure for all the order of the expansion ( |A7| ) we obtain a perturbative expression for the 
Green's function G: 



g^!{u;)\gf,{^-Jl2-qu^o)\gr{^) 



g^Uuj ~ Jl2-pujo)^gll{uj - 2J/2 -pwo - qujo)^gf'{uj - J/2-pujq) 



qO, 



-grH + ---, (AlO) 



which can be resummed in the compact form (we drop now the site indices): 



(All) 



or in the matricial form: 



5n,mjG''\Lo~nuja~ J/2). 



(A12) 



In particular the terms in square brackets in Eq. ( [A10| ) corres pond to the iteration up to the second order respectively 
of GP'"(w) and of G°^{uj - J/2 - pivo) according to Eq. (^Uj). 

Note that the contribution of the hopping processes [second term in Eq. (A12)] is diagonal in the phonon space. 
This is due to the property (A6) which relies on the locality of the electron- phonon interaction. 

From the explicit solution of the atomic problem: 



where 



we end up finally with following self-consistent equation for the Green's function in the multiphonon space: 



Lo — nujQ — -G {lu — nujQ — J/2) 



(A13) 
(A14) 

(A15) 



The solution of Eq. ( A15 ) reduces to theriaversion problem of a tridiagonal matrixEl'E l. T he diagonal elements 
G"" can be expressed as continued fractionslij obtaining Eqs. (§)-(0) and (|9|)-(pTl). Eq. (A15) looks similar to Eq. 
(34) of Ref. |l^ with the important difference that Gq^ depends now on the exchange energy J. The same result can 
be obtained for G using diagrammatic techniques as in Ref. |l^. 



APPENDIX B: ADIABATIC LIMIT 



In this Appendix we will solve the problem of one electron moving in an infinite coordination static lattice. Here 
we follow the derivation of the adiabatic limit done in the Holstein model in Ref. ^ and we use the same notations. 
The adiabatic limit is achieved as Af — > cxd keeping k — Mluq constant. The coupling constant of Hamiltonian Eq. 
(^ is given in terms of g' hy g = g' /\/2MixiQ. The polaron energy — —g'^/ojQ = —g /2k is then a well defined 
quantity in the adiabatic limit. Minimizing the ground state energy of the Holstein model with respect to the lattice 
deformation Xi around a given site i we have 



9' 



X. = ^(..). (Bl) 
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Therefore charge locahzation around a given site means also a locahzation of lattice deformations. 

The infinite coordination limit together with Eq. (Bl) implies that, for a single hole, only one site is appreciably 
distorted. Around a localization center (site 0) the nearest neighbor deformation is of the order 0(l/z), the next 
nearest neighbor deformation (0(l/z^)) and so on, so that the total charge can be spread over several shells of 
neighbors even in the z — *■ cx) limit but the deformations around a localization center vanishes in the z -t^oo limit. 

From the equation of motion we can derive a hole propagator for a given set of lattice deformations .E3 The main 
simplification of the d ~f oo limit is then that the elastic energy is solely determined by the 0-site deformation, for it 
depends on Xf. Consequently we have two kinds of local propagators: one which describes the motion of the electron 
from site back to site and which depends upon the deformation: 



and a second propagator which enters in Shopl"-^) [Eq- 
on exchange J: 

G(c.) = — 



Goo (c^) = -—^ ; (B2) 

which does not depend on lattice deformation but depends 



i^Gicu - J/2) 



(B3) 



It is worth to note that Eqs. (B2,B3) can be obtained within t — J ~ Vi model of ref.lii with on-site energy vi — g'Xo 
and zero neighbors energy. 

The lowest energy pole of Gq.o gives the electronic energy E^i for a given 0-site deformations Xq- 



E,i + g'Xo - -Re 



G 



0. 



Eq. (B4) also defines Xq as a function of E^. By defining propeply scaled deformation Xq 
Ed = te, Etot — ietot and exploiting the continued fraction structurdlZI of Eq. (B2) we have 



/ 



1/4 



J/2 



1/4 



V 



1/4 



(B4) 

g'u/k and energies 



(B5) 



e - 3 J/2 - . . . / 

By adding the elastic contribution wc have the total energy which has to be minimized with respect e: 

etot(e) = Aw2(e) - e. 



(B6) 



The total energy minimization can be carried out explicitly in th e st rong coupling limit i.e. when J/t ^ 1 or when 
A 3> 1. In these limits the continued fraction appearing in Eq. ( |B5| ) can be neglected giving a linear dependence 
for u (u — — e/2A). The minimization of Eq. ( ]B6| ) gives e = — 2A and etot = —A. This limit corresponds to a s mal l 
lattice/magnetic polaron regime. In this case the deformation "saturates" the charge deformation relation of Eq. (Bl) 
and the hole is perfectly localized on a given site. 

Another interesting case is the J/t — > limit. An analytical calculation can be done in this limit following the lines 
of Rcf. We have in this case a solution with vanishing deformation which gives the lowest energy for A < 0.844. 
In this case the Green function is the same of Ref. |l^ and consists of a semicircular band of localized states. And a a 
solution with non-zero deformation which gives the lowest energy for A > 0.844. In this case a pole emerges out of 
the band at low energies. 

The transition at A^, — 0.844 is found to be discontinuous. It is important to notice that even if it is possible to follow 
the formal steps of Ref. ^ to recover these solutions the physical interpretation of case with vanishing deformation is 
quite different. In particular we may understand this solution as describing a localized large lattice/magnetic polaron 
in the limit of infinitely large polaronic radius in contrast to the case of the pure Holstein model where in this case the 
motion of the electron is coherent through the lattice.t3 Instead solution associated to a non vanishing deformation 
has the same character in both models i.e. it describes a localized small polaron. 

In the general case the minimization of Eq. (B6) can be easily carried out numerically. Derivatives of the ground 
state energy with respect to g' and J gives respectively the hole-phonon and the exchange (mean number of spin 
defects) contributions to the total energy. The hole kinetic energy being obtained by subtraction. These derivatives 
of the ground state energy are discontinuous at the transition found for J/t — a,t \c — 0.844. The discontinuous 
large to small lattice/magnetic polaron transition exists up to J/t ~ 0.132. For larger magnetic couplings a smooth 
crossover takes place. 
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FIG. 1. Spectral density A{u!) for different values of ex- 
change energy: J/t = 0,0.2, 1.0, oo and A = 0.5, uio/t = 0.5. 
We use a finite broadening (0.02t) is used. 
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FIG. 3. Multiphonon processes in the X-J/t space for 
ujo/t ~ 0.5. From the left to the right the lines corre- 
spond to P(n = 1) > P{n = 0), P{n = 2) > P{n = 1), 
P{n = 3) > P{n — 2), etc.. Dashed lines indicate the antia- 
diabatic limit, respectively — 1,2,3, etc.. 




FIG. 4. Lattice polaron transition as determined by 
Eq. ( p^ ) in the X-uo/t space for different exchange couplings: 
J/t = 0.0 (solid line), J/t = 0.5 (dashed line), J/t = 1.0 
(dotted line). 
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FIG. 5. Moan number of spin defects A^s.d. as function of 
the exchange coupling J/t for different adiabatic parameters: 
Loo/t = 0.0,0.5,2.0 and A = (sohd hues), A = 0.5 (dashed 
hnes), A = 1.0 (dotted hnes), A = 2.0 (dot-dashed Unes). 
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FIG. 6. Effective hopping amplitude t* /t as function of the 
exchange coupling and of the electron-phonon coupling A. 
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FIG. 7. Polaron phase diagram of the Holstein-t-J model 
for different values of the adiabatic ratio wo/t- Solid 
lines mark the lattice polaron formation, dashed lines the 
large/small spin polaron transition. Inset: pictorial sketch 
of the generic phase diagram. 
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FIG. 8. Adiabatic phase diagram obtained at uio/t = 0. 
The solid line marks the large/small lattice polaron discontin- 
uous transition, the dashed line the large/small spin polaron 
transition. The sharp large/small lattice polaron transition 
disappears at Jc/t ~ 0.132 (marked by the filled circle) where 
it becomes a continuous crossover. 
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